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Abstract 

We present a combined QCD analysis of recent data produced by the HI and 
ZEUS collaborations on the diffractive and leading-proton deep inelastic positron- 
proton scattering structure functions, F^ and F^^, respectively. It is shown 
that the QCD framework for semi-inclusive processes, based on fracture functions, 
allows a unified treatment of both diffractive and leading-proton processes, and 
offers a precise perturbative QCD description for them, alternative to those that 
rely on model-dependent assumptions, such as Regge factorization, or other non- 
perturbative approaches. 
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Introduction 



In recent years considerable attention has been paid to deep inelastic 
lepton-proton scattering processes with either an identified leading-proton 
in the final state |], |2|, |, ||, or of a diffractive nature, presumably dominated 
by the single dissociative process 7*p — > Xp || [| . The two kinds of events 
differ considerably in the kinematical regions where they are produced and 
also in the non-perturbative mechanisms or models used to describe them, 
such as meson exchange in the first case, and 'pomeron+reggeon' exchange 
in the second. Consequently, their description is usually done within com- 
pletely unrelated frameworks and in terms of different structure functions, the 
leading-proton and the diffractive deep inelastic scattering positron-proton 
structure functions, and _F 2 LP< ' 3 ' ) , respectively. 

However, if the experiment under consideration is able to identify the 
final-state proton in the diffractive processes, as in the ZEUS measurements, 
or if there is confidence in the dominance of the single dissociative process 
7*p — > Xp, both kinds of events can be thought as of semi-inclusive na- 
ture, with identical final-state particles produced in the target fragmenta- 
tion region. From this point of view, the perturbative QCD framework for 
their description must be identical, with only one factorized observable, even 
though the specific models for their non-perturbative features are completely 
different. 

In perturbative QCD, the most appropriate description, for semi-inclusive 
DIS events in which the identified final-state hadron is produced in the target 
fragmentation region, is the one that includes fracture functions 0. These 
functions can be understood essentially as parton densities in an already frag- 
mented target, and there are different reasons for which they are necessary 
ingredients. In the first place, the more familiar description of semi-inclusive 
events, in terms of parton densities and fragmentation functions, refers es- 
sentially to only a limited set of current fragmentation processes; in leading 
order, final-state hadrons are taken into account only if they are produced 
in the backward direction. Going to higher orders, a breakdown of hard 
factorization is found unless certain kinematical regions are subtracted from 
the cross section in an arbitrary and process-dependent way, usually done by 
imposing cuts on the final-state variable z |§. 
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Fracture functions solve all the above-mentioned problems, introducing 
naturally the leading-order description for target fragmentation events in 
the forward direction, and also a generalization to higher orders, allowing a 
consistent factorization of collinear divergences ||. At the same time, the 
formalism provides not only the motivation for establishing a connection be- 
tween the kinematical regions of these kinds of processes, but also rigorous 
predictions about the perturbative QCD behaviour of the cross-sections. A 
priori this behaviour is not the same as that of the ordinary structure func- 
tions due to appearence of an inhomogeneous term in the Altarelli-Parisi 
equations for fracture functions |7|, [HJ . 

As for structure functions in totally inclusive deep inelastic scattering, 
QCD does not predict the shape of fracture functions unless it is known at 
a given initial scale. This non-perturbative information has to be obtained 
from the experiment, and, eventually, can be parametrized finding inspiration 
in non-perturbative models, as for ordinary structure functions. Up to now, 
the formalism of fracture functions has been successfully applied to describe 
forward neutron data from the ZEUS experiment fly]] , using a model as non- 
perturbative input. 

In this paper we present results from a QCD global analysis of recent semi- 
inclusive deep inelastic scattering data produced by HI using the framework 
of fracture functions. It is shown that this approach unifies the descrip- 
tion of diffractive and leading-proton phenomena and allows a perturbative 
QCD description without the usual assumptions about approximate Regge 
factorization. 

As a result of this analysis, we also present a parametrization for the 
fracture function that characterizes the underlying semi-inclusive process at 
a given initial scale, as extracted from HI data. The scale dependence of 
these data is found to be in remarkable agreement with the one predicted 
for the fracture function, driven by the corresponding evolution equations. 
The parametrization is also used to compute other observables measured by 
ZEUS, not included in the fit, finding also an outstanding agreement with 
the data. 

In the following section we show the relation between leading-proton and 
diffractive structure functions, introducing the concept of fracture functions 
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in this context. Then, motivated by known non-perturbative models, we 
propose a parametric form for fracture functions at an initial scale, and fix 
their parameters by means of a leading order (LO) global QCD fit to HI 
data. Finally we compute, with the resulting parametrization, semi-inclusive 
observables measured by ZEUS and present our conclusions. 



Definitions 



It is customary to define the leading-proton structure function 
from the corresponding triple-differential deep inelastic scattering cross sec- 
tion: 



d 3 a LP 4na 



i-y + ^K P(3) M 2 ,0, (i) 



dxdQ 2 d£ xQ A 
with the usual kinematical variables 

W " q X ~2P-q 4 " q-P P ~ 2q-(P-P')~ C [ ' 

Here q, P, and P' are the momenta of the virtual boson, of the incident 
proton, and of the final state proton, respectively, and y = Q 2 /(xs). Both 
the cross section and the structure function defined in this way, include an 
implicit integration over a given range for the small transverse momentum 
of the final-state proton. 

Even though the processes accounted for in eq. (1) are explicitly of a 
semi-inclusive nature, the formulation based on the leading-proton structure 
function is used instead of the usual approach for semi-inclusive deep inelastic 
scattering 



d 3 a?.„, r . r . Ana 1 



l-2/4F£ e . 2 *(^ 2 ) x ^^ 2 ) ( 3 ) 



dxdQ 2 dz xQ A \ 2 

in terms of parton distributions and fragmentation functions, because the last 
one only takes into account hadrons produced in the current fragmentation 
region and thus not contributing to the forward leading hadron observables. 

These problems, and those related to factorization of collinear singular- 
ities at higher orders, are overcome, however, if the complete perturbative 
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framework for semi- inclusive processes is taken into account, for which the 
LO expression for the production of very forward hadrons is [] 



d 3 aL r „^ 47ra 2 



,2 N 



T """ ! " ll-y+ y -)j:e*xM?^(x,z,Qi), (4) 



dxdQ 2 dz xQ 4 I " 2 



i 



where Mf /p (x, % z, Q 2 ) is the fracture function that accounts for target frag- 
mentation processes. 

Fracture functions can be thought of in terms of the elements of any 
lowest-order picture for hadronization, for example as the product of a flux 
of exchanged 'reggeons' times their structure functions, or more formally 
as an ingredient of the perturbative QCD treatment: the non-perturbative 
parton distributions of a proton fragmented into a proton. The latter choice 
can be generalized to higher orders, allowing a consistent analysis of the 
scale dependence of the semi-inclusive cross section. This scale dependence 
is driven by inhomogeneous Altarelli-Parisi equations 



dM^x,z,Q 2 ) a s (Q 2 



dlogQ 2 2vr 

a s (Q 2 ) f x /( x+z ) dy 



Jx/(i-z) y 3 \y J 



(5) 



2vr 




where Pij(y) and Piji(y) are the regularized and real splitting functions, re- 
flecting the fact that the evolution is not only driven by the emission of 
collinear partons from those found in the target (homogeneous evolution), 
but also by the fragmentation of partons radiated from the one struck by the 
virtual probe (inhomogeneous evolution). 

Defining the equivalent to F2 for fracture functions, i.e. 

M p 2 /p (x,z,Q 2 ) =*£efMf/ p (*,*,Q 2 ), (6) 

i 

and taking into account the shift from z to £, the relation between this func- 



tion and the leading-proton structure function is quite apparent. In ref. [Ill 



1 For the sake of simplicity we define the variable z = Ep> /Ep as the ratio between the 
energies of the final-state proton and the proton beam in the centre of mass of the virtual 
photon-proton system. For very forward protons then, 1 — £ ~ z. 
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the use of fracture functions for the description of leading-hadrons produced 
in the target fragmentation region has already been discussed in relation to 
the analysis of very forward neutrons observed by the ZEUS collaboration. 
There it was shown that, neglecting contributions beyond LO coming from 
the current fragmentation region, what is usually defined as the leading- 
neutron structure function f£ N ® is just the fracture function contribution 
to the semi-inclusive cross section. 

Analogously to eq. (1), the differential cross section for diffractive deep 
inelastic scattering is usually written in terms of the diffractive structure 
function 



where x p = C,, and the variable (3 is used instead of x. As in the previous 
case, an integration over the small transverse momentum of the final-state 
proton is implied. In this context it is customary to use the variable t = 
(P — P') 2 instead of the transverse momentum. When this last cross section 
is dominated by the single dissociative process YP Xp, implying that 
there is a proton in the final state, the contributions to it are again given 
by the same fracture function in (4), even though in a completely different 
kinematical region^]. The diffractive region is given by small values of x p 
(xjp < 0.1), whereas leading-proton data are associated with larger values 
of x p ( x p > 0.1). The available data also differ considerably in the ranges 
covered by the variable (3, 0.04 < (3 < 0.9 in the case of diffractive structure 
functions, and very small values ((3 < 0.03) for leading-proton data. 

Different kinematical regions correspond to different behaviours and also 
to different underlying models. The leading-proton structure function has 
been measured by the HI collaboration || [|], and has also been compared 
with predictions of different mechanisms, such as meson exchange and soft 
colour interaction models for example, implemented in event generator pro- 

2 The integration over the measured range of the variable t can straightforwardly be 
circumvented by defining 'generalized' p"2| fracture functions with an explicit dependence 
on that variable and obeying homogeneous evolution equations, analogously to the diffrac- 
tive structure function F^ 4 \ However, in the present analysis we restrict ourselves to 
ordinary fracture functions. 
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grammes [O, ITJI], none of which reproduced integrally the main features of 



the data ||. 

The standard interpretation for the diffractive cross section is given in 
terms of 'pomeron' exchange. In this framework, different model estimates, 
and even QCD- inspired parametrizations of the 'pomeron' content, have been 
proposed |l5[ 16, 17]. More recently a formidable set of data has been pre- 



sented by the HI collaboration, and as in the leading-proton case, the com- 
parison between model predictions and data has been found to be rather 
poor unless a large number of additional elements is included P . 

In the language of fracture functions, both the leading-proton and the 
diffractive regimes are complementary features of a more general semi-inclusive 
process. The approach, then, suggests and provides a bridge between the two 
regimes, which is particularly appropriate and even necessary, at least in the 
kinematical region where neither the 'pomeron' nor the 'reggeon' exchange 
picture alone are expected to describe the data, say x p ~ 0.1. 



Parametrization 



Having introduced the arguments for a combined analysis of both leading- 
proton and diffractive cross sections using the framework of fracture func- 
tions, we proceed in this section to their phenomenological determination by 
means of a LO QCD global analysis of the available dataQ. 

In order to obtain a parametrization for the proton-to-proton fracture 
function Mf /p (/5, Q\, x p ) at a given initial scale Qq, we select in the first 
place, a relatively simple functional dependence in the variables (3 and x p , 
but with enough flexibility as to reproduce the data accurately. If one were 
only interested in the diffractive regime, the natural choice would be a simple 
'pomeronic' flux in x p times an ordinary parton distribution in (3, namely 

M p 2 lp {^Qlx p ) ~C p x a p x N 8 a °(l-P) b °, (8) 

where the label s stands for singlet contribution, assumed to be the dominant 

3 Since the data were obtained under the assumption of a negligible longitudinal diffrac- 
tive structure function, we keep the analysis at leading order for consistency. 
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one. Analogously, for the leading-proton regime the natural choice would be 
almost the same, but with a standard meson or 'reggeon' flux or even better, 
something combining their effects [TS| : 



M% /p (P,Ql,x lp )~C LP x a p ^a-x 1P ) f5LP x N s (3 a °{l-l3) b *. (9) 

These kinds of parametrizations give relatively good initial approxima- 
tions to the description of the corresponding data sets; however their survival 
seems unlikely in a more precise analysis. For example the 'pomeron' flux 
factorization hypothesis (independence with respect to /3 and Q 2 ) is a rather 
strong assumption, not even true for the kinematical range of the available 
data. This has been been shown in recent analyses |J and in fact is apparent 
from fig. 1 where there is a clear x a behaviour for large values of /3, with 
a ~ —0.25, behaviour that is then suppressed for intermediate values, and 
finally is changed to something like a positive power of x F at low (3. 

In order to take into account small departures from the initial approx- 
imations, and also combine the two behaviours in such a way that for low 
Xjp (diffractive regime) the 'pomeron' picture dominates, while for low f3 and 
large x p the meson or 'reggeon' exchange picture emerges, we propose a 
modified flux such that the light-quark singlet component (M p ' p = 3M P ' P = 

3M^ P = 3Mf/ p ) of the fracture function is parametrized as 

xM^(P,Ql,x p ) = iV s /^(l-/3) 6 * x 

{c P (3 x y + c LP (i - py™ (1 + a LP (l - x p )^)} , (io) 

and similarly for gluons with the corresponding parameters N g , a g and b g . 

Notice that this is only done for the sake of convenience; even though at 
the initial scale Ql the parametrization implies some sort of flux factorization, 
beyond the initial scale, the evolution equations drive the fracture function as 
a whole making the usual discrimination between 'fluxes' and 'parton den- 
sities of the exchanged object' somewhat ambiguous. This is particularly 
relevant in this case, where none of the 'pomeron' and the 'reggeon' contri- 
butions are known, making it then possible to shift part of the pomeron-like 
contribution into the other during the fitting procedure by simply adjusting 
the parameters conveniently. 
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Even at the initial scale, the flux ambiguity is reflected in the fact, that for 
parametrization purposes, the three normalizations, N s , C p , and Clp, are 
not independent, and consequently, any of them can be set to unity (we take 
C p = 1). We do not assume any momentum sum rule for the distributions 
in this analysis, since its significance is not clear even when considering the 
process in terms of the exchange of 'pomerons' and 'reggeons'. 

Regarding the evolution, we set Ql = 2.5 GeV 2 , take A QC d = 0.232 GeV 2 , 
and we choose to work in a scheme with a variable number of flavours, where 
charm and bottom distributions are radiatively generated from their corre- 
sponding thresholds (taken as their masses). The main effect of the heavy 
quarks in the fit is to reduce the light-quark singlet component, effect that is 
particularly sizeable for charm quarks due to the smaller mass (almost coin- 
cident with the initial scale), and the influence of the large gluon distribution 
in the evolution. 

In the analysis of both HI data sets 0, ||, which were obtained within dif- 
ferent ranges of the variable t, we have assumed the universal validity of the 
exponential behaviour in t measured by the ZEUS collaboration || and have 
rescaled the data to a common range, given by that of the diffractive data of 
HI (— MpX 2 p /(l — x F ) = —t min < — t < 1 GeV 2 ). The parametrization ob- 
tained through the fitting procedure corresponds to that particular range of 
the variable t; it has to be rescaled for the computation of any other observ- 
able if different kinematical cuts are applied. This is particularly important 
for the leading-proton HI data, since a small transverse momentum cut was 
implemented px < 200 MeV, selecting only data corresponding to \t\ <C 1 
GeV 2 . Such a small cut leads to a large rescaling factor (close to 3). We 
have also taken into account an additional error of 5% for the leading-proton 
data, because the uncertainties in the rescaling process. 

Phenomenological parametrizations like the one in eq. (10) lead to si- 
multaneous QCD fits for both the diffractive and the leading-proton cross 
sections obtained by HI with very good precision. In order to ensure that 
the evolution code is working properly in the region of large (3, we imposed a 
lower constraint on the (1 — (3) exponents for both quark and gluon distribu- 
tions, b s , b g > 0.1. Furthermore, in the expectation of very hard distributions 
we saturate this constraint and fix a s = (leaving only 8 free parameters). 
Doing this we obtain a global fit with Xtota//d-o.f.=1.09 (x 2 = 292.23, data 
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points = 274), which we designate by Fit A. 

Similar parametrizations, but with softer gluons, yield slightly higher 
values, for example in Fit B, b g is set to 0.7 (whereas b s remains 0.1) finding 
Xi 0to ;/d.o.f.=1.13. The values found for the parameters of both sets at the 
initial scale Ql are given in table 1. 



Set 




Clp 


I3lp 


Ilp 


dLP 


N s 




a g 


A 


-1.260 


14.395 


32.901 


2.627 


12.320 


0.041 


0.354 


0.450 


B 


-1.257 


12.556 


32.412 


2.338 


11.412 


0.047 


0.694 


0.648 



TABLE 1: Parameters for Ql = 2.5 GeV 2 . 

As is shown by the solid lines in figs. 1 and 2, the accuracy of the fit is re- 
markably good in the case of HI diffractive data (x 2 H i D /data = 215.63/226), 
and also good for the leading-proton data. The dashed line in fig. 1 shows 
the contribution coming from the x F F term in eq. (10), which could be 
interpreted as the 'pomeronic component' of the fracture function. This con- 
tribution is clearly dominant for large (3, but fails to do so in the low-/3 and 
large-x^ region where the contribution from the second term in eq. (10) 
changes the x p slope of the distributions. 

In fig. 3, we show the light-quark singlet and gluon fracture densities at 
the initial scale Ql and at an intermediate value of Q 2 = 10 GeV 2 for two 
characteristic values of x p . The first one, x p = 0.005, which corresponds 
to well inside the diffractive regime, shows a rather hard behaviour, whereas 
the second one, x p = 0.2, which belongs to the leading-proton regime, is 
much softer. It is worth noticing that, at large x p , the large-/? behaviour 
of the distributions is not well constrained because leading-proton data are 
only available in the very small-/? range. 

As can be seen from the plots, the differences in the initial distributions 
between fits A and B tend to be washed out by the evolution and are mostly 
sizeable for the gluon at /3 > 0.5. As also shown in fig. 3, gluons carry much 
more impulse than quarks, specially in the case of the diffractive regime and 
at low scales. The evolution damps down this gluon dominance, suggesting a 
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large fraction of valence-like gluons in the 'pomeron', which is not so apparent 
for 'reggeons'. Putting both mechanisms together in the physical fracture 
function (fig. 4), the ratio between the 'momenta' carried by gluons and 
quarksQ depends strongly on the value of x p , which makes one or the other 
dominant. Even though gluons clearly dominate, parametrizations with only 
gluons at the initial scale, or even at lower values (Qq = 1 GeV 2 ), and thus 
having only quarks of radiative origin, lead to global fits of considerably 
poorer quality (x 2 /d.o.f. > 1.5) 

The main conclusion that can be drawn about the gluon density from the 
analysis is that the distribution is important at large 0, at variance with the 
one for inclusive structure functions, but its behaviour cannot be precisely 
determined yet from the available data, in particular, the exponent of the 
(1 — 0) factor in the parametrization. 



Scale dependence 



As it has been said, the fracture function approach leads to very definite 
predictions about the scale dependence of the semi-inclusive cross sections. 
The rigorous factorization of the cross sections achieved within this formalism 
allows a precise perturbative QCD analysis of the scale dependence of the 
data, as is usually done for ordinary structure functions. 

In figs. 5 and 6 we compare HI diffractive and leading-proton data, at 
fixed values of x p and 0, as a function of Q 2 , with the evolved fracture 
function. The scale dependence induced by the evolution equations (5) is 
perfectly consistent with the data. For very large values of (0 ~ 0.9) 
there seems to be some evidence of higher-twist contributions, as has been 



predicted within some perturbative models [19 



It is also worth noticing that within the measured range, the scale de- 
pendence is dominated by the homogeneous evolution. The effects of the 
inhomogeneous term become important only for low values of xl where frag- 



4 This ratio is defined as the one between the integrals of x M g and x M q in the range of 
(3 [0, 1]. In this way, when cither of the exchange pictures dominates, the ratio coincides 
with the usual definition for the fraction of the momentum carried by partons in cither 
the 'pomeron' or the 'reggeon'. 
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mentation functions are larger. These effects are, however, beyond the range 
of present data, as can be seen in fig. 7. There, the solid lines represent 
the scale dependence induced in the fracture function by ordinary evolution 
equations (homogeneous term in eq. (5)), whereas the dashed lines take into 
account both contributions. In order to compute the inhomogeneous contri- 
bution we have used GRV94 parton distributions | 2Cj , and the fragmentation 
functions of ref. [21]. 



ZEUS measurements 



In order to check the distributions obtained in the previous section, and 
also our assumption about the t-dependence, in the following we compare 
results obtained using our best parametrization (Fit A) with data sets pre- 
sented by the ZEUS collaboration, which had not been included in the fit. 

First, we compare ZEUS diffractive measurements || with the outcome of 
the parametrization after the appropriate evolution to the mean scale value 
of the data Q 2 = 8 GeV 2 , and the already mentioned rescaling in t. Figure 
8 shows the agreement found between the data and the fracture function 
estimate for different bins in (5. 

It is important to stress that at variance with HI diffractive data, the 
ZEUS data are obtained by identifying the final-state proton, i.e. explicitly 
collecting semi-inclusive events. The agreement found gives us confidence in 
the assumption about the dominance of single dissociative processes 7*p — > 
Xp in HI diffractive data, and then the applicability of a semi-inclusive 
framework. 

Going to leading-proton data, the ZEUS collaboration has measured the 
fraction of DIS events with a leading-proton in the final state giving special 
attention to the dependence in the variable t [|j]. At least two different sets 
of data have been presented, corresponding to different ranges in t. After 
the adequate rescaling of the parametrization for the fracture function at an 
average value of Q 2 = 10 GeV 2 , and using GRV94 pi parton distributions 



for the estimate of the total number of events [I1J , we obtain a remarkable 
agreement with the data in the common x^-range, as shown in figs. 9a and 
9b. The Q 2 dependence is found to almost cancel in the ratio between the 
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corresponding fracture and structure functions. 

Again, the general framework, and particularly the assumption about 
the t dependence, seems to be nicely confirmed. Notice also the fact that 
the parametrization interpolates fairly well the diffractive and leading-proton 
behaviour between HI kinematical regions (thick lines) where the data used 
in the fit belong. 

The resulting interpolation is radically different from those obtained, for 
example, adjusting diffractive data with a purely 'pomeronic' flux, indicat- 
ing the relevance of other contributions in the region x p ~ 0.05-0.1. These 
effects are of extraordinary importance when addressing the comparison be- 
tween diffractive DIS and pp collision data as a mean of testing factorization 

0- 

Conclusions 

We have shown that an approach based on fracture functions motivates, 
from a formal point of view, and also phenomenologically allows a unified de- 
scription of both leading-proton and diffractive deep inelastic cross sections. 
A simple parametric form for this function gives a very accurate description 
of the data available at present providing a smooth interpolation between 
the distinctive behaviours of the two regimes, also in accordance with ZEUS 
data. The analysis also hints at some non-perturbative features, such as a 
strong-gluon dominance in the 'pomeronic' component with a characteristic 
valence-like behaviour but with still large uncertainties in the large-/? region. 
Finally, our results verify that the scale dependence of the data agrees with 
the one predicted by the fracture function formalism. 

A code containing both sets of parametrizations can be requested from 
the authors by e-mail at daniel.deflorian@cern.ch or sassot@df.uba.ar 
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Figure 1: HI diffractive data against the outcome of the fracture function 
parametrization (solid lines) and its pomeron-like component (dashed lines). 
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Figure 2: HI leading-proton data against the outcome of the fracture function 
parametrization (solid lines). 
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Figure 3: Fracture function densities: a) and b) light-quark singlet, 
c) and d) gluons. 
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Figure 5: Scale dependence of HI diffractive data and the one obtained 
evolving the Fit A parametrization (Fit B renders undistinguishable results). 
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Figure 6: Scale dependence of HI leading-proton data and the one obtained 
evolving the Fit A parametrization (Fit B renders undistinguishable results). 
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Figure 7: Homogeneous and complete evolution of the proton-to-proton 
fracture function (solid and dashed lines, respectively). 
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Figure 8: ZEUS diffractive data, against the expectation coming from the 
fracture function parametrization (Fit A). 
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Figure 9: The fraction of DIS events with a leading-proton for different 
ranges of t as measured by ZEUS, against the expectation coming from 
the fracture function parametrization (Fit A). 
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